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We consider statistical mechanics models of continuous spins in a 
disordered environment. These models have a natural interpretation 
as effective interface models. It is well known that without disorder 
there are no interface Gibbs measures in infinite volume in dimen- 
sion d — 2, while there are "gradient Gibbs measures" describing an 
infinite-volume distribution for the increments of the field, as was 
shown by Funaki and Spohn. 

In the present paper we show that adding a disorder term prohibits 
the existence of such gradient Gibbs measures for general interaction 
potentials in d = 2. This nonexistence result generalizes the simple 
case of Gaussian fields where it follows from an explicit computation. 

In d = 3 where random gradient Gibbs measures are expected to 
exist, our method provides a lower bound of the order of the inverse 
of the distance on the decay of correlations of Gibbs expectations 
w.r.t. the distribution of the random environment. 

1. Introduction. 

1.1. The setup. Our model is given in terms of the formal infinite- volume 
Hamiltonian 

(1) H[r]]{ip) = - ipj)-J2viVi- 

Here the fields (unbounded continuous spins) € M represent height vari- 
ables of a random surface at the site z € Z"^. Such a model is motivated as 
an effective model for the study of phase boundaries at a mesoscopic level 
in statistical mechanics. The disorder configuration -q = {r]i)i^^d denotes an 
arbitrary fixed configuration of external fields, modeling a "quenched" (or 
frozen) random environment. 
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For background and various earlier results about both continuous and 
discrete interface models without disorder, see [2, 5, 8, 12, 13] and refer- 
ences therein; for results about discrete interface models in the presence of 
disorder, see [3, 4]. 

Assumptions. The pair potential V{t) is assumed to be even, V{t) = 
V{—t), and continuously differentiable. We require that V grows faster than 
linearly to infinity, that is, lim|f||oo jip^r = oo for some positive e. p{-) is the 

transition kernel of a simple random walk on Z*^, assumed to be symmetric 
and of finite range. We further demand for simplicity that the random fields 
rji be i.i.d. under the distribution P, symmetric and with finite nonzero 
second moment. We denote the expectation w.r.t. P by the symbol E. 

Vector fields, fields and gradient fields. We call the set of sites Z"^ with 
oriented edges between i,j whenever p{i — j) > the graph of the random 
walk. We call a vector field (or vector field configuration) a map from the 
set of oriented edges ij of this graph to M such that Vij = —yj,i. Every 
field configuration U = (f/i)jeZ'* gives rise to a vector field configuration 
U' = U'{U) by Ulj = U{i) — U{j). In this case we call U' the gradient field 
of U. There does not need to exist such a function for a general vector field 
V. Its existence is equivalent to the loop condition Vij + Vj^k + yk,i + Vi,! = 
along plaquettes i,j,k,l. We denote the set of all field configurations in 
infinite volume by and the set of all gradient field configurations in infinite 
volume by 0'. 

Gibbs measures and gradient Gibbs measures. The quenched finite- volume 
Gibbs measures (or local specification) corresponding to the Hamiltonian (1) 
in a finite volume A C Z'^, a boundary condition ip and a fixed disorder con- 
figuration in A are given by the standard expression 
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where Z'^[ri] denotes the normahzation constant that turns the last expres- 
sion into a probabihty measure. It is a simple matter to see that the growth 
condition on V guarantees the finiteness of the integrals appearing in (2) for 
all arbitrarily fixed choices of rj. 

We note that the Hamiltonian H[ri] changes only by a configuration- 
independent constant under the joint shift (px*-^ ^x + coi all height variables 
with the same c. This holds true for any fixed configuration ij. Hence, finite- 
volume Gibbs measures transform under a shift of the boundary condition 
by a shift of the integration variables. Using this invariance under height 
shifts we can lift the finite-volume measures to measures on gradient config- 
urations, defining the gradient finite-volume Gibbs measures (gradient local 
specification). These are the probability kernels from 0' to fi' given by 

(3) J (/i')i'N(V)G'(/) := / fil[v]{dip)G{ip'{^)), 

where (p is any field configuration whose gradient field is (f' . 

Finally, we call fx' a gradient Gibbs measure if it satisfies the DLR equa- 
tion, that is, 

(4) J J {p')f{d^')f{p') = J ii'{dp^')f{^'). 

Here we dropped the i] from our notation. 

Gradient Gibbs measures have the advantage that they may exist, even 
in situations where a proper Gibbs measure does not. When this happens, 
it means that the interface is locally smooth, although at large scales it has 
too large fluctuations to stay around a given height. For background about 
these objects, we refer to [8, 9, 12]. 

Translation- covariant {gradient) Gibbs measures. Denote by T^if = 
{tpi-r)iizid the shift of field configurations on the lattice by a vector r S U^. 
A measurable map rj ^ ^[rj] is called a translation- covariant random Gibbs 
measure if ^[r/] is a Gibbs measure for almost any and if it behaves ap- 
propriately under lattice shifts, that is, / ^[Trri\{dip)F{ip) = J ii[i]]{d{p)F{Trip) 
for all translation vectors r. This means that the functional dependence of 
the Gibbs measure on the underlying disordered environment is the same 
at every point in space. This notion goes back to [1] and generalizes the 
notion of a translation-invariant Gibbs measure to the set-up of disordered 
systems. 

Naturally, a measurable map rj is called a translation- covariant 

random gradient Gibbs measure if / fi'[Trr]]{dip')F{ip') = J fi'[r]]{dip')F{Tr(p') 
for all r. 
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1.2. Main results. A main question to be asked in interface models is 
whether the fluctuations of an interface that is restricted to a finite volume 
will remain bounded when the volume tends to infinity, so that there is 
an infinite-volume Gibbs measure (or gradient Gibbs measure) describing a 
localized interface. This question is well understood in translation-invariant 
continuous-height models, and it is the purpose of this note to discuss such 
models in a random environment. 

Let us start by discussing the nonrandom model in the physically inter- 
esting dimension d = 2. We start with the Gaussian model where V{t) = 
Gaussian models are simple because explicit computations can be done. 
These show that the Gibbs measure ^^{d'^) does not exist in infinite volume, 
but the gradient field (gradient Gibbs measure) does exist in infinite vol- 
ume. Both statements may easily be derived by looking at the properties 
of (differences of) the matrix elements of (/ — -Pa)"^ which appears here as 
a covariance matrix, where Pa is the transition operator, restricted to A, 
given in terms of the random walk kernel p. Equivalently, one may say that 
the infinite-volume measure exists conditioned on the fact that one of the 
variables (pi is pinned at the value zero. Funaki and Spohn showed more 
generally that for convex potentials V there are tilted gradient Gibbs mea- 
sures and, conversely, a gradient Gibbs measure is uniquely determined by 
the tilt [8, 9, 12]. 

For (very) nonconvex V , new phenomena are appearing: There may be 
a first-order phase transition in the temperature where the structure of the 
interface (at zero tilt) changes, as shown by Biskup and Kotecky [2]. This 
phenomenon is related to the phase transition seen in rotator models with 
very nonlinear potentials exhibited in [6, 7]. The basic mechanism is an 
energy-entropy transition such as was first proven for the Potts model for a 
sufficiently large number of spin values [10]. 

What can we say for the random model? In [11] the authors showed a 
quenched deterministic lower bound on the fluctuations in the anharmonic 
model in a finite box of the order square root of the sidelength, uniformly in 
the disorder. In particular, this implies that there will not be any disordered 
infinite- volume Gibbs measures in d = 2. This latter statement is not sur- 
prising since there is already nonexistence of the unpinned interface without 
disorder. But what will happen to the gradient Gibbs measure that is known 
to exist without disorder, once we allow for a disordered environment? 

Gaussian results and predictions. Let us first look in some detail at the 
special case of a Gaussian gradient measure where V{t) = specializing 
to nearest neighbor interactions. Then, for any fixed configuration r/A, the 
finite-volume gradient Gibbs measure with zero boundary condition is a 
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Gaussian measure with expected value 

y&K 

(5) 

where T^^ = (-Aa)-^^ - (-Aa)tI. 

The matrix coy f,{ipij]Lpim) = {-^Kji} - {-^kXli - (-'^^Ui + (-^a)^^ 
is the same as in the model without disorder; its infinite- volume limit exists 
in any dimension d>2, by a simple computation. 

What about the infinite- volume limit of the mean value (5), as a function 
of the disorder? We note that X^[r]] := {X-^It]]) itself a random 

vector field with covariance 

(6) CAiij, Im) = n^l) ^ Tl^,,yTit^y. 

j/GA 

Moreover, it is also a gradient vector field, since the loop condition carries 
over from ip' , by linearity of the Gibbs expectation. 
For its variance, we have, in particular, 

(7) CK{ij,ij)=nril)Y.iTiy?- 

y€A 

In two dimensions the infinite-volume limit of this expression does not exist 
since r log r)^ dr ~ log N, when the sidelength N of the box diverges to 
infinity. In dimension d > 2, we have r'^~^{-^r~^'^~'^^)^ dr r~^'^~^^ dr, 
so the fluctuations stay bounded. 

In particular, this explicit computation shows that in the Gaussian model 
there cannot be infinite- volume random gradient Gibbs measures in d = 2. 
Indeed, already the local — short-distance — fluctuations are roughening up 
the interface. 

It is the main result of this paper to show that this result persists also for 
an anharmonic potential where explicit computations are not possible. 

Theorem 1.1 (Nonexistence in d = 2). Suppose d = 2. Then there does 
not exist a translation- covariant random gradient Gibbs measure fj.'[r]](dip') 
that satisfies the integrability condition 



(8) E 
for sites i ^ j. 



< oo 



Let us consider dimension d = 3 where translation-covariant infinite- volume 
gradient measures are believed to exist. The next result shows that, if they 
do, they must have slow decay of correlations w.r.t. the random environment. 
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Theorem 1.2 (Slow decay of correlations in d = 3). Suppose that d = 3 
and that fi'[r]]{dip') is a random gradient Gibbs measure. Put 

(9) C{ij,kl) :=E^J ^x'Md^')V'{^',^) J fi[rj]{dip')V'{^',, 

for sites i ^ j and k ^l. 
Then 

(10) lim sup ' /' = C30 
'•T~i,j,fc,h|i-fc|>r r y^+^i 

for all e > 0. 



Remark. Note that in d = 3 for the case of the quadratic nearest neigh- 
bor potential V{t) = y a translation-covariant random gradient Gibbs mea- 
sure exists by explicit computation. It is the Gaussian field whose covariance 
and mean are given by the infinite- volume limits of (6) and (5). An explicit 
computation with this Gibbs measure also shows that the power of the bound 
given in the last theorem is optimal; see Section 2.3. 



The method of proof of both theorems relies on a surface- volume compar- 
ison. It is inspired by the Aizenman-Wehr method (devised in [1] and used 
on discrete interfaces in [4]), but different in several aspects. As opposed to 
the Aizenman-Wehr Ansatz, where free energies are considered, we derive a 
discrete divergence equation r] = VX where the external random field (the 
disorder) acts as a source and the vector field X is provided by the expecta- 
tion of V'{(pi — ipj) with respect to the hypothesized gradient Gibbs measure. 
Exploiting a discrete version of Stokes' theorem and a volume versus surface 
comparison between terms, we arrive at a contradiction. 



2. Proof of Theorems 1.1 and 1.2. 



2.1. Proof of Theorem 1.1. The method of proof is to argue from the ex- 
istence of a translation-covariant gradient Gibbs measure to a contradiction 
in d = 2. To do this, we start with the following definition. 

Definition 2.1. Let be an infinite- volume gradient Gibbs measure, 
in a random or nonrandom model. Then we call the vector field X on the 
graph of the random walk given by 

(11) X,J:=J^,'{d^'){V'{^[^)) 

the associated vector field. 
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This indeed defines a vector field because, by symmetry, V'{x) = —V'{—x) 
and, hence, Xij = —Xji. Of course, the same definition can be made in finite 
volume, but for our proof, we will work immediately in infinite volume. 

We note that a gradient Gibbs measure in the Gaussian model provides 
even a gradient vector field, since V'{x) = ex is a linear function, and so 
the loop condition carries over from (p'^j to Xij. For general nonquadratic 
potentials V, the vector fields X^j will not be gradient fields. This explains 
why, in the anharmonic model, we must work with vector fields (functions 
on the edges). 

Proposition 2.2. Let fi'[r]]{dip') be a random infinite-volume gradient 
Gibbs measure and Xij[ri] := J fi'[r]]{d(p'){V'{ip'^j)) the associated vector field 
on the random walk graph for i ^ j. Then 



for all i € IJ^ . This equation will be called the ^^divergence equation^\ 

Remark. This equation is a discrete version of the equation 77 = V • X 
which holds in continuous space W^. 

Comparing with Maxwell's equations, X plays the role of the electric field, 
and r] plays the role of the electric charge. In the Gaussian case, we know 
that X is a gradient-field, that is, curl-free. In the case of general potentials 
V it will not be. Note that the kernel of this equation, that is, the vector 
fields X which are the solutions of the equation = V • X, consists of all 
divergence-free vector fields, which are plenty. In particular, the divergence 
equation does not allow to determine X uniquely in terms of 

Proof of Proposition 2.2. We look at the one-site local specification 
at the site i. Take a single-site integral over ipi appearing in the partition 
function and use partial integration to write 



(12) 




(13) 
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Now, integrate over the remaining (pk for in a finite volume A to see that 

(14) 

j 

Integrating this equation over the boundary condition (p' w.r.t. fJ-'[i]] and 
using the DLR equation for the gradient measure implies the proposition. 
□ 

Summing the divergence equation over i in a finite volume A, we note 
that the contributions of the edges that are contained in A vanish, due to 
the property of X being a vector field. Hence, we arrive at the following 
corollary. 

Corollary 2.3 (Integral form of divergence equation). 

(15) J2V^= E PU-i)X^M■ 

ieA 

Note that the sum of boundary terms on the r.h.s. plays the role of a 
surface integral in the Stokes equation. 

Let us now specialize to d = 2 and prove Theorem 1.1. We choose A = 
{— L, —L + 1, . . . , L}^ and normalize by j^. 

We remark that the sum over boundary bonds ij decomposes into the 4 
sides of a square, whose bonds will be denoted by Bl{u)^ v = 1,2,3,4. 

So we have 

i£A ieAjigA'^ 

(16) 

By the ergodic theorem, each of the four sums converges to its expected 
value 

(17) liml Pii-j)X^M = '^ E PUrnXoAv]) 

in whatever sense (e.g., almost surely or in L^), where 61,62,63 = —61,64 = 
—62 are unit vectors on Z'^. 
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By the CLT, the l.h.s. of equation (16) does not converge almost surely 
to a constant; indeed, j^J^ieAVi converges only in distribution to a Gauss 
distribution with a nonzero variance. This is a contradiction to convergence 
to a delta-distribution, which proves Theorem 1.1. 

2.2. Proof of Theorem 1.2. Take the square of the integral form of the 
divergence equation and take its expectation w.r.t. the measure P: 

(18) nVo)\M= E E p{j-i)p{k-l)Q,,ki. 

ieAjGA=A:eA,/GA': 

Assume the bound 

(19) Qj-fc«<Const(l + |i-A:p)-'?. 

Let us take for A a ball w.r.t. the Euclidean metric of radius L. Then, for 
large L, the l.h.s. behaves like a constant times L^. 

The large-L asymptotics of the r.h.s. of (18) is provided by the large-L 
behavior of the double integral 



(20) / dX{x) dXivXl + lx-ylT", 

where LS^ denotes the two-dimensional sphere with radius L and A is the 
Lebesgue measure on the sphere. By rotation-invariance, this equals 

(21) cL^ f dX{y){l + \y-Le\^)-'i = cL^ ( dX{z){l + L^\z - efy^ , 

where e is the North Pole of . Using polar coordinates cos 9 = s for the 
point z, where s = 1 is corresponding to the point e, we have \z — ep = 
2(1 — s). This gives 



(22) 



/ dX{z){l + L'^\z - ef)-i 

/I 97r9l~9 
ds{l + ^^2(1 - s)y = — (1 + 2L^ 
-1 ^ - Q 



By (18), we thus have under the assumption of (19) that < cL'^~^'^, for 
large L, which can only be true for q < ^- This implies that Cij^kl can not 
decay faster than the inverse distance between i and k. 

2.3. Sharpness of polynomial decay in d = 2>. 

Proposition 2.4. Suppose d='i. Let ^'[rj\{d(p') he a random gradient 

2 

Gihhs measure for the Gaussian model V{x) = ^ with i.i.d. rji with finite 
second moment. 
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Then 

(23) lim sup ' ^ ^' <oo. 

'•T~ij,fc,i:|i-fc|>r 

Remark. This shows that the fluctuation lower bound is sharp in d = 3. 

Proof of Proposition 2.4. Let us choose i = —Re and k = Re, where 
e is a unit coordinate vector so 2R is the distance between i and k. To 
estimate the large- i2 asymptotics of the decay of the covariance in infinite 
volume which is given by 

(24) C^3 {ij, Im) = E{r,l) J] TflyT^^y, 

s/ga 

we are led to consider the large- i? asymptotics of the following integral: 

1 1 



(25) I{R) := J d'y. 



l + \y-Re\^l + \y + Re\^' 

Now the proposition follows by an explicit computation. 

Indeed, using cylindrical coordinates, this integral can be rewritten as 



re 

I{R) = 2tt 

Jo 



dz I dff- 



1 + (Z - i?)2 + f2 1 + (z + i?)2 + f2 



^ ~4RJo ~s^Jo "^'[l + iz-Ry + s' l + {z + R)^ + s 

= Ir'^^^- 

Here we have used the substitution s = and 

f^dz^ l + (z + Rf f°°du^ R-'^ + {u + lf 

27 J{R):= —log— 7 —log 7^2-7 TX2- 

Jo z l + {z-R)'' Jo u R ^ + {u-iy 

Finally, the function J{R) is increasing in R and has the finite limit 

(28) lim J{R) log = 8 ^ = vr^. 

^ ^ moc ^ ' Jo u ^(n-l)2 Jo x2-l 

(See an integral table or check with Mathematica.) This shows the proposi- 
tion. □ 
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